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In order to study the processes creating intermediate and heavy nuclei in massive stars it is 
necessary to provide neutron capture cross sections and reaction rates close to stability and for 
moderately unstable neutron-rich nuclei. Furthermore, one has to know the efficiency of neutron- 
releasing reactions in the main evolutionary phases of a massive star. We present simulations of the 
nucleosynthesis in a 15 and 25 Mq star, for the first time followed completely from main sequence 
hydrogen burning until the type II supernova explosion including all nuclides up to Bi. Theoreti- 
cal reaction rates were calculated with the NON-SMOKER code, providing a complete library of 
Hauser-Feshbach cross sections and rates for nuclear and astrophysical applications. Experimental 
rates at stability were taken from different sources. The impact of uncertainties in the rates on 
nucleosynthesis are illustrated by two examples, the reactions 62 Ni(n,7) 63 Ni and 22 Ne(a,n) 25 Mg. 

I. INTRODUCTION 

Stars above ~ 10M Q are responsible for producing most of the oxygen and heavier elements found in nature. 
Numerous studies have been devoted to the evolution of such stars and their nucleosynthetic yields, e.g., by Refs. J|] 
(in the following referred to as WW95) or However, our knowledge of both the input data and the physical processes 
affecting the evolution of these stars has improved dramatically in recent years. Thus, it became worthwhile to improve 
on and considerably extend the previous investigations of pre- and post-collapse evolution and nucleosynthesis. We 
present the first calculation to determine, self-consistently, the complete synthesis of all stable nuclides in any model 
for a massive star. 

We employed a nuclear reaction network of unprecedented size in full stellar evolution calculations. The network 
used by WW95, large in its days, was limited to 200 nuclides and extended only up to germanium. Studies using 
reaction networks of over 5000 nuclides have been carried out for single zones or regions of stars with simplified 
networks, especially to obtain the r-process, e.g., but "kilo-nuclide" studies of nucleosynthesis in complete stellar 

models (typically of 1000 zones each) have been hitherto lacking. With the continuous improvement in understanding 
the input physics as well as in computer technology over the last several years it became worthwhile and feasible to 
re-evaluate the nucleosynthesis in massive stars using most recent physics and numerics. 

In the following discussion we focus on the nuclear physics inputs and their relevance to the production of interme- 
diate and heavy nuclei. Due to the increasing Coulomb barrier reactions involving neutrons play an important role, 
either by directly producing the weak s-process component and in a possible n-process (i.e. an intermediate process 
between s- and r-processes), or by providing the seed nuclei from which the 7-process produces proton-rich isotopes 
via photodisintegration processes. 



II. THE MODEL 

Our calculations were performed using the stellar evolution code KEPLER [jj]j6| with several modifications regarding 
astrophysical inputs and the numerical treatment. The latter includes introduction of a novel adaptive reaction 
network for nucleosynthesis considering all relevant nuclides from hydrogen thru polonium at any given time; this 
amounts to about 2350 nuclides during the explosion. The largest extension of the network is shown in Fig. [j]). For 
details see It is evident that there are many isotopes far off stability included for which experimental information 
is scarce. This is even more true in more exotic scenarios such as the r-process which involve isotopes close to the 
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FIG. 1. Maximum size of reaction network during explosive burning. The line of stability is indicated. 

dripline. (In our calculations we do not follow the proposed r-process in the ^-wind emanating from the proto-neutron 
star shortly after the collapse of the Fe core.) 

As in previous studies |Q] the explosion was simulated using the piston model. However, the piston procedure was 
slightly modified. We also included the ^-process changing the light element abundances during the explosive phase 
with the same cross sections as in WW95. We updated the opacity tables (OPAL95) |§| and considered mass loss due 
to stellar winds fl(i[| . 

III. THEORETICAL AND EXPERIMENTAL NUCLEAR INPUT 

Here, we only briefly summarize the most important nuclear inputs. For details on the complete rate set, see Ref. 

A. Statistical Model Rates 

As can be seen in Fig. [l], a large part of our network is comprised of unstable nuclei for which no experimental 
information is available. Therefore we have to heavily rely on theoretical predictions. It can be shown that 
the majority of those reactions proceeds via the compound nucleus reaction mechanism and can be described in the 
framework of Hauser-Feshbach theory. As the basis for the creation of our reaction rate set we used statistical model 
calculations obtained with the NON-SMOKER code |l2|JT^ 1. A library of theoretical reaction rates calculated with 
this code and fitted to an analytical function - ready to be incorporated into stellar model codes - was published 
recently JljJ. It includes about 30000 n-, p-, and a- induced reaction rates for all possible targets from proton- to 
neutron-dripline and between Ne and Bi, thus being the most extensive published library of theoretical reaction rates 



80 — 



60 — 



N 

40 — 



20 — 



2 




3 



to date. It also offers rate sets for a number of mass models which are suited for different purposes. For the network 
described here we utilized the rates based on the FRDM set. For other ap plications, nuclear cross se ctions for the 
whole mass range are published separately [[l4| and are also available online ||http: / /data. nucastro.org | . 

Of special interest are a capture reactions on isospin self-conjugated (N — Z) target nuclei. The cross sections of 
such reactions are suppressed by isospin effects and require special treatment in theoretical models. This is described 
separately in Ref. 



B. Weak Rates 

An important change of the weak interaction rates for A < 40 is brought about by recent work flilf . For the effects 
of those new rates on the stellar structure results, see Ref. |0. The theoretical f}~ and (3 + rates of Ref. were 
also considered. As a special case, we implemented a temperature-dependent 180m Ta decay ||l9|| . Updated neutrino 
loss rates were also used |po|| but we did not follow the ^-process for nuclides with Z or N larger than 40. 



C. Experimental Rates 

The above described theoretical rates formed the backbone of the huge set of reaction rates that was required to 
follow nucleosynthesis in our extended reaction network. Whereever possible, the theoretical rates were substituted 
by or supplemented with experimentally determined reaction and decay rates. 

The experimental /3~ , (3 + , and a decay rates and branching ratios given in Ref. |^lj were implemented, including 
decays of the first excited states. 



Neutron capture rates were inferred from a compilation of 30 keV Maxwellian Averaged Cross Sections \22\ by 
renormalizing the NON-SMOKER rates to the given values. 

The most important reaction is the reaction 12 C (a,j) 16 O. Rs value at 300 keV has still not be determined satis- 
factorily in experiments (this is called the single most important problem in experimental nuclear astrophysics) but 
a variation of this rate impacts the C/O ratio obtained from He burning and thus the stellar structure and the total 
subsequent evolution and nucleosynthesis of the star. The rate previously used in WW95 was the one of Ref. |23| ] 
multiplied by 1.7. Here, the more recent evaluation of J24|] was used which yielded ,9(300) = 146 keV barn. The 
temperature dependence was chosen as suggested by p5| |. 

As a standard rate for 22 Ne(a,n) 25 Mg we used the rate of Ref. (2(| which is close to the lower limit of the NACRE 
rate. This is also the rate used for the results shown in Fig. 0. The uncertainties in this rate are discussed below in 



Sec. [VB. 



IV. SELECTED IMPORTANT NUCLEAR REACTIONS 

The results from our 15 M Q model are shown in Fig. ||. The obtained abundances are not only sensitive to details 
in the stellar physics (see e.g. J2?]]) but also to nuclear physics. Nuclei beyond the Fe peak can in general only be 
produced, directly or indirectly, via neutron capture reactions which underlines the fact that it is important to have 
a sound knowledge of the nuclear processes induced by neutrons with energies of a few tens of keV. 

In this section two important reactions involving neutrons are discussed in more detail. Both give rise to considerable 
uncertainties in the final isotopic yields obtained in our model. They are not the only reactions of special interest but 
due to the limited space we want to use these to illustrate the difficulties one has to deal with when implementing 
nuclear reaction rates in nucleosynthesis calculations. The first reaction - 62 Ni(n,7) 63 Ni - is important to model 
a specific detail in the abundance distribution whereas the second reaction - 22 Ne(a,n) 25 Mg - is the main neutron 
source and therefore impacts the synthesis of all intermediate and heavy nuclides. 



A. The Reaction 62 Ni(n,7) 63 Ni 

Examination of the production factors of the Ni isotopes in Fig. || clearly shows a pronounced overproduction of 
62 Ni in our calculation. This especially catches the eye because all other nuclides are produced in amounts nicely 
agreeing with the band defined via 16 O. A closer look at the problem reveals a peculiar situation for the reaction 
62 Ni(n,7) 63 Ni which destroys 62 Ni. For the results shown here, the reaction rate was taken from Ref. [Q. Interestingly, 
by comparison with the previous version |2£| of that reference it is found that the rate has changed by about a factor 
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FIG. 3. Production factors in a 25 M© star with NACRE recommended rate for 22 Ne(a,n) 25 Mg. 

of three between these two compilations. The old value is higher by about a factor of three which leads to increased 
destruction, sufficient to bring the 62 Ni production down to a value comparable to the other nuclides in that region 
in our calculation of the 15 M Q star. 

Searching for the reason for the change in the 62 Ni(n,7) 63 Ni rate (cross section), it turns out that the new value 
is not based on a more recent experiment but rather on a different extrapolation of the cross section measured only 
at thermal neutron energies. The astrophysically relevant neutron energies are in the range 20 < E n < 40 keV and 
the Maxwellian Averaged Cross Section at 30 keV is tabulated in Refs. [p2pq ]. Since there is only a measurement of 
62 Ni(n t h,7) 63 Ni using thermal neutrons n t h, an extrapolation to 30 keV is needed. The older value uses a standard 
1/v extrapolation, typical for direct capture of an s-wave neutron. However, it was realized p9fl that the 14.5 barn 
thermal neutron capture cross section of 62 Ni is unusually high as compared to its neighboring isotopes, those showing 
thermal cross sections which are lower by factors of 4 — 6. This lead to the idea that the thermal 62 Ni capture cross 
section is enhanced due to a subthreshold resonance in addition to the direct capture contribution |^9|| . More specific, 
the excess value was attributed to the high energy tail of a bound s-wave resonance at —0.077 keV. Since the resonance 
contribution decays faster than 1/v to higher energies, at 30 keV only the contribution of the direct term survives 
which is then lower by a factor of 3 than the term assumed in the older rate extrapolation. Unfortunately, neither the 
actual strength of the subthreshold resonance nor the neutron capture cross section at 30 keV have been determined 
directly by an experiment. 

It would be audacious to infer from the discrepancies in our calculations that the actual rate of 62 Ni(n,7) 63 Ni is lower 
by about a factor of 3 than the new value given in Ref. and that the previously quoted value would be preferrable. 
However, in the light of the experimental situation it is justified to direct the attention of the experimentalists to that 
problem and to point out the importance of neutron capture measurements in the astrophysically relevant energy 
range instead of at thermal energies. 
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FIG. 4. Production factors in a 25 Mq star with upper limit of the NACRE rate for 22 Ne(a,n) 25 Mg. 

B. The Reaction 22 Ne(a,n) 25 Mg 

This reaction is the main neutron source in the type II supernova model. It provides the neutrons to build up 
intermediate and heavy s-process nuclei in the secondary s-process component originating from type II SN. These 
nuclei are the seed nuclei for the 7-process which produces proton-rich, stable nuclei (p-nuclei) by photodisintegration 
(see e.g. |l] j30| , |3l|| ). The produced p-nuclei can be nicely seen as the enhanced proton-rich isotopes sticking out from 
the floor of weak s-process abundances for A > 90 in Fig. |^. Thus, a change in this rate impacts a large number of 
nuclei, basically all of the nuclides beyond Fe. 

It has also been speculated that a largely enhanced 22 Ne(a,n) 25 Mg rate could provide enough neutrons to make an 
r-process (e.g. Refs. ||3^-|3^|). Despite of the exponential dependence of the rate on stellar temperature, this would 
require extreme temperature conditions, not found in our models. At most, we find slightly increased temperatures 
during supernova shock front passage through the He shell, leading to a modest increase in the neutron flux over the 
standard s-process flux. We call this an n-process which leads to a slight redistribution of abundances but produces 
no r-process features and no significant total yields. 

When trying to implement a proper rate for 22 Ne(a,n) 25 Mg one finds that the experimental errors are still con- 
siderably large. For example, comparing the results when using the recommended value for that reaction from the 
NACRE compilation [^H| and the upper limit quoted in the same compilation, we find that the production factors are 
higher by factors of 10 — 30 for all nuclides beyond Ga (Figs. ||, |^). Recently, it was claimed that such a high value 
could cure the underproduction of the p-process isotopes in the Mo-Ru region which has been found in all stellar 
models j3(|. Indeed, also those p-nuclei are produced more strongly, however, in a consistent calculation like ours one 
clearly sees that such an enhancement is ruled out by the fact that other nuclides are equally enhanced. This not only 
limits the possible maximum value of the 22 Ne(a,n) 25 Mg rate but also shows that another solution has to be found 
for the production of the light p- nuclides in amounts comparable to the more heavy p- nuclides. 

Again, this shows the importance of measuring the cross sections in the relevant energy range which is difficult in 
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the case of 22 Ne(a,n) 25 Mg. Because of its importance in heavy element nucleosynthesis it is sometimes called the 
second-most important reaction in nuclear astrophysics, after the well-known reaction 12 C(a,7) 16 0. 



V. SUMMARY 



The first calculation to self-consistently determine the complete synthesis of all stable nuclides for massive stars 
was performed with a reaction network of unprecedented size. Examples for the results for a 15 M Q and a 25 M Q 
star were shown above. 

We find a strong sensitivity of the production of intermediate and heavy nuclei on the 22 Ne(a,n) 25 Mg rate which 
still remains to be measured to satisfactory accuracy in the relevant energy range. 

As another example we pointed out that even neutron capture reactions on stable isotopes, which are supposed 
to be well known, are still a source of large uncertainties. A number of those reactions have only been measured 
at thermal neutron energies and were then extrapolated to 30 keV. This is always problematic because the energy 
dependence of the cross section can well be different from the usually assumed 1/v behavior, due to p-wave capture 
and possible contributions of the tails of resonances, both for resonances at higher excitation energies and subthreshold 
resonances. It has to be strongly suggested that the (n,7) databases should be searched for such extrapolated cases 
and that experimental efforts are undertaken to measure across a wider range of energies. This would greatly improve 
the reliability of the nuclear input and the credibility of the results of nucleosynthesis calculations. 
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